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The method of equations of motion for spin operators in the case of 0(3) Heisenberg ferromagnet 
is systematically analyzed starting from the effective Lagrangian. It is shown that the random phase 
approximation and the Callen approximation can be understood in terms of perturbation theory for 
type B magnons. Also, the second order approximation of Kondo and Yamaji for one dimensional 
ferromagnet is reduced to the perturbation theory for type A magnons. An emphasis is put on 
the physical picture, i.e. on magnon-magnon interactions and symmetries of the Heisenberg model. 

Calculations demonstrate that all three approximations differ in manner in which the magnon- 
magnon interactions arising from the Wess-Zumino term are treated, from where specific features 
and limitations of each of them can be deduced. 
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I. INTRODUCTION 

The study of magnon-magnon interactions is a prob¬ 
lem with long history [1] whose complexity initiated de¬ 
velopment of numerous analytical tools. Its importance 
goes far beyond pure academic interest since with the ad¬ 
vent of new materials, especially those connected to the 
high-temperature superconductors [2-5], the question of 
magnon-magnon interactions and their influence on ther¬ 
modynamic properties of ordered magnets became an ur¬ 
gent one to solve. This amounts not only to finding effi¬ 
cient calculation techniques, but also understanding how 
or why certain approximations work and what are their 
domains of applicability. With this in mind, it is use¬ 
ful to study simpler models [6]. In this case, they reveal 
some features generic to the wide class of spin systems 
otherwise concealed in more complicated ones contain¬ 
ing second and third neighbor interaction, various types 
of frustration and anisotropies, random impurities etc. 
Much like 2D Ising model, due to several well established 
and exact results, the (quantum) 0(3) Heisenberg ferro¬ 
magnet has come to be a prototypical model for testing 
diverse theoretical and numerical methods. 

Nearly all theories addressing the problem of magnon- 
magnon interactions are spin-operator oriented in the 
sense that all calculations are based on the operators 
used to define the Heisenberg model (See Fig. 1). The 
most straightforward way to incorporate the interac¬ 
tion effects is to express the localized spins in terms of 
bosonic/fermionic operators directly in the Heisenberg 
Hamiltonian or to use the coherent-state path integral 
(See e.g. [7, 8]). The original spin Hamiltonian is after¬ 
wards interpreted as describing the system of interacting 
bosons/fermions, to which standard perturbation theory, 
or suitable mean-field approximations (MFA), may be 
applied (See [I, 7-9] and references therein). Although 
techniques build upon boson/fermion representations of 
spin operators provide an important insight into the be¬ 
havior of magnetic systems, they suffer from a universal 
disadvantage: The su(2) algebra of spin operators and 
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the dynamics of the spin system are fully satisfied only 
with exact boson/fermion Hamiltonian and correspond¬ 
ing Hilbert space. Accordingly, any approximation in the 
boson/fermion Hamiltonian, such as approximate expres¬ 
sion for localized spins or some mean-field approximation, 
destroys the spin nature of operators, often in 

an uncontrolled manner. 

A different way to deal with magnon-magnon inter¬ 
actions is based on the equations of motion (EOM) for 
spin operators [10]. Aside from the fact that the EOM 
for spin operators are practically unsolvable for general 
Heisenberg-type Hamiltonian, the analysis through EOM 
gets more involved giving the fact that spin operators 
present basic dynamical variables and charge densities 
at the same time [11]. Despite all that, the method of 
equations of motion, especially in temperature-Green’s 
function (TGF) disguise (See [12] for a recent review and 
original references), gained quite a popularity. There 
are several reasons for this. First, its predictions for 
thermodynamic properties of Heisenberg-type magnets 
agree with Monte Garlo simulations [12] and with exper¬ 
imental results (See e.g. [13-17] and references therein). 
Equally important is the flexibility of the method mak¬ 
ing it easily adjustable to systems with complex [18-22] 
or low-dimensional lattices [23-30], second and third- 
neighbor interaction or anisotropies [13, 15, 31-35], which 
is of great importance when studying real compounds. 
Also, the TGF method is recognized as useful in theo¬ 
ries of diluted magnetic systems [36, 37], nuclear spin or¬ 
der in quantum wires [38], multiferroics models [39, 40] 
and even in theories of itinerant electron systems where 
Heisenberg Hamiltonian appears as an intermediate ef¬ 
fective model [22, 41, 42]. The approximations made di¬ 
rectly in the equations of motion, known as the decou¬ 
pling schemes (DS) in the TGF jargon [12], enable one 
to solve linearized system of EOM-s thereby determin¬ 
ing the magnon spectrum, correlation functions etc. By 
suitably choosing the parameters of linearization [43, 44] , 
more or less satisfactory results may be obtained. Even 
though TGF method is extensively exploited in contem¬ 
porary research (See the papers quoted earlier in this 
paragraph), a thorough discussion on EOM and magnon- 
magnon interactions seems to be lacking. This gap in the 
literature comes from the fact that the nonlinear Hamil- 
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FIG. 1. Different perspectives on dynamical degrees of freedom and magnon-magnon interactions in 0(3) ferromagnets: 

a) In spin-operator-based approaches, interactions are viewed as a consequence of nontrivial algebra and Hilbert space of spin 
operators coupled through exchange integral J [One-dimensional lattice of localized S = 1/2 spins is displayed for clarity] 

b) In EFT, magnons are considered as a set of (local) coordinates on the sphere = 0(3)/0(2) and magnon-magnon 
interactions are dictated by the underlying geometry of 5^, the non-Abelian character of 0(3) and the existence of conserved 
charges in ferromagnetic ground state. 


Ionian describing magnon-magnon interactions need not 
be introduced in EOM approach, which makes the com¬ 
parison with usual boson/fermion representation theories 
rather difhcult. 

The brief sketch of boson/fermion representation and 
EOM for spin operators from two preceding paragraphs 
deliberately points on methodical differences in han¬ 
dling magnon-magnon interactions, as they appear to be 
among frequently used analytical tools and thus cover a 
significant part of approximate treatments of spin sys¬ 
tems. Equally important, however, is their common fea¬ 
ture: Being based on spin operators and not on the true 
degrees of freedom (i.e. magnons), these approximations 
tend to obscure the physics behind them. This is es¬ 
pecially true for linearized EOM, since calculations are 
conducted without direct reference to magnon-magnon 
interaction operators. 

An alternative view on magnons in a system of local¬ 
ized spins, based upon symmetry considerations, is put 
forward in [45]. As an extension of the phenomenologi¬ 
cal Lagrangian method [46-52] to the nonrelativistic sys¬ 
tems, the outcome of [45] is the effective Lagrangian for 
the 3D 0(3) Heisenberg ferromagnet constructed from 
the magnon fields and Tr'^(x) transforming nonlin- 

early under 0(3) rotation group (The subject of non¬ 
linear realizations and effective Lagrangians is reviewed 
in [53], while texts with special emphasis on nonrela¬ 
tivistic systems include [54-58]). The main advantage 
in viewing magnons as Goldstone bosons, accompany¬ 
ing 0(3) —0(2) symmetry breakdown, is that their 
dynamics is completely separated from the su(2) alge¬ 
bra of spin operators entering the Heisenberg Hamil¬ 
tonian. Magnons are considered as local coordinates 
on 0(3)/0(2) = and all information about magnon- 
magnon interactions are encoded in the nonlinear La¬ 
grangian (See Eig. 1). The perturbative calculations may 
be conducted in a systematic manner, free of the upper 
mentioned ambiguities characterizing spin-operator ori¬ 
ented approaches. The general program of [45] is suc¬ 
cessfully applied on the calculation of the free energy up 
to two [59, 60] and three loops [61] and is subsequently 
extended to low dimensional ferromagnets [62-65] . If one 


is concerned only with low-temperature series of the free 
energy and related quantities, which are strongly con¬ 
trolled by internal symmetry of the Heisenberg model, 
the continuum field-theoretic methods of [45, 56, 59-65] 
will suffice. However, the comparison with experimen¬ 
tal data and other spin-operator-based methods often 
requires for the discrete symmetry of the lattice to be in¬ 
cluded also. This is the case, for example, with magnon 
energies whose dispersion deviates from simple law 
due to the lattice anisotropies. The symmetry of the lat¬ 
tice may be incorporated in the effective field description 
by the method of [66] (See also [56]), resulting in an ef¬ 
fective theory which describes ferromagnet with the help 
of lattice magnon fields. This means that the internal 
0(3) symmetry is implemented through the usual non¬ 
linear realization (equivalently, symmetry may be real¬ 
ized linearly on the magnon fields, but with an additional 
nonlinear constraint), while the lattice symmetries enter 
via adequately chosen regulator. At the same time, the 
use of spatial lattice singles out canonical formalism in 
constructing quantum field-theoretic description of a fer¬ 
romagnet. Beside the fact that the Hamiltonian lattice 
theory and the functional Lagrangian method with di¬ 
mensional regularization yield the same results, as it was 
shown in [66] , few other important issues favor canonical 
over functional integral formalism in this particular case. 
Eirst, with the lattice magnon Hamiltonian at hand, a 
comparison with standard methods based on boson rep¬ 
resentations of the spin operators [1, 7, 8, 67, 68] becomes 
a much easier task. Second, EOM method (i.e. two-time 
temperature Green’s function theory) which is be sys¬ 
tematically reexamined in present paper, is intimately 
related to the spin operator Hamiltonian through the 
definitions of thermodynamic averages. Also, the pro¬ 
cedure of quantization itself is insightful regarding how 
magnon-magnon interactions may appear in an effective 
field theory. And finally, interacting Hamiltonian for fer¬ 
romagnetic magnons, rather than the general form of in¬ 
teracting Lagrangian, reveals certain simplifications in 
loop calculations [66] if a Lagrangian contains the Wess- 
Zumino term. 

The lattice Hamiltonian theory of magnon fields pro- 
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vides us with an efficient tool for studying effects of 
magnon-magnon interactions in 0(3) ferromagnets, and 
likewise, the interactions induced by linearizations in 
EOM approach. This is an extremely important point 
since, according to the current understanding [ 12 ], the 
quality of DS can be judged only a posteriori. Also, 
it is difficult to distinguish between the true effects of 
magnon-magnon interactions induced by linearizations of 
EOM from pure artifacts of TGF formalism. All this 
can lead to incorrect interpretation of results obtained 
using EOM for spin operators. Previous remarks bring 
us to the main purpose of the present article. By ex¬ 
plicitly identifying types of magnon-magnon interactions 
that actually constitute some of the well known DS and 
connecting them with internal symmetry of the 0(3) fer- 
romagnet, we put EOM approach in broader context of 
interacting spin waves and effective Lagrangians thereby 
connecting it with recent developments in these fields. 
At the same time we clarify some misinterpretations of 
EOM method in the case of 0(3) ferromagnets. Quite 
generally, the effects of magnon-magnon interactions can 
be inferred either from low-temperature expansions of 
the free energy and related thermodynamic quantities 
or from renormalized magnon energies, as correspond¬ 
ing results in linear spin wave (LSW) approximation are 
well known [ 66 ]. Since the heart of EOM method lies 
in determination of the renormalized magnon energies, 
we shall concentrate on finding appropriate effective La- 
grangian that yields, through the perturbation theory, 
desired magnon self-energy. This will allow us to de¬ 
termine the true nature of some well known lineariza¬ 
tions from the physical point of view, i.e. relying only 
on magnon-magnon interactions and internal symmetry 
of the Heisenberg model. 

The perturbation theory for lattice magnon fields 
shares many features with the standard effective field 
theory. The most important distinction concerns calcu¬ 
lation details, as the spatial lattice replaces continuum 
field theoretic setting. Section II summarizes main defi¬ 
nitions and relations needed for loop calculations which 
are presented in Sec. Ill and IV. Specifically, it is shown 
in Sec. Ill that the random phase approximation (RPA) 
and the Callen approximation (CA) can be understood 
as interaction theories for type B magnons (see [57, 69] 
for A/B classification of Goldstone fields). Finally, Sec. 
IV describes the second-order approximation of Kondo 
and Yamaji for one dimensional ferromagnet as a pertur¬ 
bation theory for type A magnons. Further justification 
of the lattice magnon field formalism comes from com¬ 
parison with Quantum Monte Carlo simulation [sections 
III F and IV C 3] which, at the same time, supports some 
of our conclusion regarding RPA, CA and KYA. 


II. LATTICE LAPLACIANS AND 
MAGNON-MAGNON INTERACTIONS 

Conventional effective field theory (EFT) combines 
derivative and loop expansions to set up perturbative cal¬ 
culations to given order in some energy scale [53-55] . The 
relevant scale for ferromagnets is determined by tempera¬ 
ture and EFT enables one to systematically calculate free 


energy in powers of T [51, 59-65]. Yet, as emphasized in 
the Introduction, some issues concerning ordered mag¬ 
nets are adequately treated by the perturbation theory 
for lattice magnon fields which preserves the full symme¬ 
try of Heisenberg Hamiltonian. This approach, adjusted 
for 0(3) ferromagnets and systematically exposed in fol¬ 
lowing sections, includes couplings of neighboring sites 
through the lattice Laplacian. For general lattice it is 
defined by [70] 




on _ Y 

5^y:[«x+A)-0(x) 


( 1 ) 


where x denotes a lattice site, D is the dimensionality of 
spatial lattice and {A} are vectors connecting each lattice 
site with its Zi nearest neighbors. The definition (1) can 
be extended in an obvious manner to next-nearest neigh¬ 
bor couplings. The plane wave solutions, corresponding 
to free magnons, are built with eigenvalues of 
2D 

exp[ifc • x] = -T-rry [1 - lD{k)\ exp[ifc • x] (2) 

l-'l 

= —k^ exp[ifc • x], lD{k) = exp[ifc • A]. 

IM 


Depending the presymplectic structure [57, 69] , Goldsone 
bosons may be classified as type A or B. In the present 
context, magnons of type A posses free energies propor¬ 
tional to k and those of type B to . As a consequence 
of derivative couplings, the loop integrals generally con¬ 
tain eigenvalues combining loop and external momenta, 
p — q^. These can be simplihed according to 

J{nq)o ^ = j (%)o P ^ + g ^ 


^2 ^2 

P Q 


, (3) 

Iq Jq 

where (np)o denotes Bose distribution for free magnons 
and 

f-f 


/ — / (27r)^' 

The focus in Sec. IV is on one dimensional ferromagnets 
described by an effective theory that includes couplings of 
lattice magnon fields spanning beyond nearest neighbors. 
To make presentation in Sec. IV more transparent, we 
shall denote with V^ 2 ) the Laplacian for second neighbor 
couplings on a chain. It is defined by 


[1 - 7(2g)] = - 2q 


so that 


2q =g 2 


1 « -2 
l--q 


(5) 


( 6 ) 


and 7 (g) is defined in (2) for D = 1. An analogue of (3) 
for V( 2 ) reads 


/ {nq)o 2 p - 2q 
J a 




2q 


— 2a^ 2p 2q 


(7) 


Equations from this section will be used throughout the 
paper for loop calculations, frequently without direct ref¬ 
erence. 
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III. PERTURBATION THEORY FOR TYPE B 
MAGNONS 

We introduce lattice magnon Hamiltonian and accom¬ 
panying diagrammatic rules for type B magnons [57, 69] 
in this section. They will serve as a basic for the sys¬ 
tematic investigation of several approximate treatments 
of magnon-magnon interactions in 0(3) ferromagnes for 
H > 3. 


A. Lattice magnon Hamiltonian 


The starting point in an effective field description of 
a ferromagnet is the classical Lagrangian [45] (See also 
[71-73]) of the unit vector field f/® := [H^, C/^, = 

[ 7 r(a;), [/^(a;)]'^ transforming linearly under 0(3) 


-b EnHU^. 


( 8 ) 


S = NS/V is the spontaneous magnetization per unit 
volume at T = OK, is a constant, to be fixed latter, 
H is the external field and a refers to the spatial coordi¬ 
nates. We also recall that the Lagrangian ( 8 ) collects all 
magnon-magnon interactions up to p^. Due to the non¬ 
vanishing charge densities in the ground state and the 
non-Abelian character of the 0(3) rotational group, the 
effective Lagrangian contains appropriate Wess-Zumino 
(WZ) term. Being linear in temporal derivatives, WZ 
term is responsible for nonrelativistic dispersion of fer¬ 
romagnetic magnons [45, 74, 75] and ( 8 ) describes the 
system of type-B Goldstone bosons [69]. In contrast to 
this well known fact, its influence in perturbation theory, 
i.e in loop corrections to the magnon self-energy and the 
ferromagnet free energy, is not fully appreciated. The 
reason for this is that WZ term is practically invisible in 
approaches based on boson representations, dominating 
the literature on magnon-magnon interactions [1, 67, 68 ], 
since it is hidden inside the commutation relations of spin 
operators. 

As outlined in the Introduction, we shall employ 
canonical formalism in constructing the quantum theory. 
We first note that the field describi ng p hysical magnons 
is complex [45, 59, 69], = \/S/ 2 [ 7 r^(a;) -|- i 7 r^(a;)], 

since the equation governing their dynamics is of a 
Schrodinger type. The standard canonical prescription 
[76] then yields the canonical momentum in the form 
n = 2i'!/)/[I + U^\, with = sJX — (2/E)'!/)t'!/>. This 
path towards consistent quantum theory is of course com¬ 
pletely legitimate, but nonlinear connection between 
and n sends important part of the magnon-magnon in¬ 
teractions from the Hamiltonian 'H('0,^^) to the canon¬ 
ical commutation relations, ['(/'(a;),n(y)] = id(a; — y). 
The quantum field theory of this sort is strongly rem¬ 
iniscent of the initial, spin-operator oriented, approach 
and is in fact something we would like to avoid. For¬ 
tunately, all complications concerning canonical com- 
mutaion relations can be circumvented by noting that 
terms with single temporal derivative may be eliminated 
from the Lagrangian ( 8 ) with the help of the equation 


of motion for U{x,t) [76] (See also [77]). In this man¬ 
ner, we obtain the free magnon Lagrangian (bilinear in 
magnon fields tt^ and tt^) and nonlinear part describ¬ 
ing their interaction [ 66 ]. Canonical quantum interac¬ 
tion theory may be constructed starting from this La¬ 
grangian. For the present purposes, we shall be needing 
terms up to and including eight magnon field operators. 
An straightforward calculation yields redesigned Hamil¬ 
tonian [ 66 ]. This Hamiltonian, however, includes only 
magnon-magnon interactions of the order p^. It is well 
known [49, 51, 52, 60, 61, 78] that consistent loop ex¬ 
pansion must include interaction terms of higher order 
in magnon momenta. Following examples from Lorentz- 
invariant chiral perturbation theory [79, 80], we include 
higher order terms in momentum, i.e. in spatial deriva¬ 
tives, by putting the Hamiltonian on the (spatial) lattice. 
We find terms with four and eight magnon field operators 

Heff = Ho + Hint, (9) 

Ho = --^vo'^'tpHx)V'^ilj{x), 

2mr) 

X 

E 

^int — 21^ + rl^ + + ilg , 

where V{) denotes the volume of the unit cell and 

H^"^ = ^vo'^7t'^{x)n{x) ■V'^7 t{x), ( 10 ) 

X 

nf^ E 7r^(a;)VV2(a:), 

X 

H^^ = ■ VV(a;), 

X 

Hs^ = -^vo E ^ ^ : ( 11 ) 

X 

We have set H = 0 for the time being. In what fol¬ 
lows, tjj will always be written to the right in expres¬ 
sions like TT ■ TT. We also note that E = S'/t;o in a lat¬ 
tice theory, with S denoting the magnitude of localized 
spins. The LSW spectrum is fully recovered if we choose 

= JS'ZijAj^/(2D) = l/(2mo) [see equation (13) 
bellow]. This simple connection between parameters of 
the effective theory and those of original Heisenbeg model 
is a consequence of the ferromagnetic ground state. 

Some comments on interacting magnon Hamiltonian 
(9)-(ll) are in order. We observe that magnon-magnon 
interactions come in two categories, distinguished by su¬ 
perscripts a and b in (lO)-(ll). While and Hg^^ 
are typical interaction terms for theories of unit vector 
fields, those containing operator 7r(a;) • V^7r(a;), namely 

and , describe the magnon-magnon interactions 
originating in WZ term. Thus, they are specific to the fer¬ 
romagnetic system and careful treatment of the magnon- 
magnon interactions arising in WZ term is mandatory, if 
the correct description of a ferromagnet is to be expected. 
We shall come back to this point several times in the text. 
It can also be seen from (9) and (10) that Ho+H^f’'^ 
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is quite similar to the Dyson-Maleev Hamiltonian [81-83]. 
Both the Dyson’s construction and the Maleev boson rep¬ 
resentation [84] yield magnon Hamiltoninan equivalent 
to i?o -I- -I- , but only the method based on the 

effective Lagrangian (9) traces 7r^(x)7r(a;) • V^7r(a:) inter¬ 
action back to the WZ term. Finally, the Hmiltonaian 
given in (9) and (lO)-(ll), being the Hamiltonian of true 
Goldstone fields, makes the weakness of magnon-magnon 
interactions at low momenta rather obvious. Thus, 1/5 
expansion, a usual starting point in theories based on 
boson representations and a priori justification for the 
weakness of magnon-magnon interactions [1, 7, 67, 68] 
is, in fact, unnecessary. The magnons interact weakly at 
low momenta irrespective of the spin magnitude 5, and 
this is what actually makes the perturbation theory well 
defined [85] even for S = 1/2. 


B. Diagramatic description of magnon-magnon 
interactions 


Before proceeding with concrete perturbation theory 
calculations, we pause to introduce a useful variant of 
Feynman diagrams [66]. It is particularly suited to the 
magnon Hamiltonian (9), but it can be easily adapted 
to other scalar field theories with derivative couplings. 
This method will allow us to systematically keep track of 
individual contributions in various magnon-magnon pro¬ 
cesses, which will turned out to be of great importance 
when comparing magnon perturbation theory with ap¬ 
proximations in the equations of motion. To account for 
the eight-magnon vertices, we slightly modify the nota¬ 
tions used in [66]. 


We shall denote with 


4(a) 


IM], and 


/m\ the vertices corresponding to interaction terms 
from (lO)-(ll). Depending on the structure of the vertex, 
lattice Laplacians may act on one, two or four magnon 
propagators. Apart from simplest, one-vertex diagrams, 
higher order corrections will always contain at least two 
lattice Laplacians acting on propagators. The calcula¬ 
tions soon become quite involved, since the number of 
propagators, as well as the number of lattice Laplacians 
acting upon them, increases. To make them feasible, 
we associate with each vertex a color so that the ver¬ 
tex and the propagators affected by its lattice Laplacian 
are drawn using that color. The rest of the lines are sim¬ 


ply drawn black. We label each propagator with D 1 
dimensional momentum k := [k,uJn]'^ and the colored 
ones also carry —k^. Here —k^ denotes the eigenvalue 
of the lattice Laplacian, defined in (2), and = 27rnT 
are the Matsubara frequencies. When colored line passes 
through the vertex, the momentum in —k^ corresponds 
to algebraic sum of the incoming and outgoing momenta. 
If more than one Laplacian acts on a certain propagator, 
the line will carry one color for each vertex. We shall 
adopt the following convention on colored propagators 
belonging to a closed loop: If the lattice Laplacian of 


4(a) 


or 



acts on a magnon propagator clos¬ 


ing a loop around it, the line will b e draw n half-colored. 
When the Laplacian comes from 4(5) or A(FK , 


the line will be drawn in full color if the Laplacian acts 
on both contracted magnon field operators. Otherwise, 
the line will also be drawn half-colored. With this con¬ 
ventions we generalize the method of [66] to interacting 
terms with more than four magnon operators. The dia¬ 
grams with colored lines and vertices will be referred to as 
the colored contractions. Later, in Sec. IV, we shall allow 
more than one Laplacian per vertex. Lastly, the magnon 
propagator and the free one-magnon energies for type B 
magnons are given by 


D{x - y,T^ - Ty) = (T {tl}{x,T^)ip'^{y,Ty)])Q (12) 

1 “ r Qiq-{x-y)-iuj„{r^-ry) 


wo(9) 


2mo ’ 


(13) 


and /3 ^ Z]n[‘^o(p)-iw„] ^ = (np)o, where (np)o denotes 
Bose distribution for free magnons. 

As an explicit example, consider a three-loop correc¬ 
tion to the magnon self-energy arising from For¬ 

mally, it is represented by 



Several colored contractions contribute to this particular 
diagram. These are 
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2mo 


'qpQ 


{nq)o{np)o{nQ)o p + q - Q - 


(17) 


since two colored loops closed around single vertex, as 
well as single colored closed loop with two colored exter¬ 
nal lines, vanish due to momentum conservation; 



This is the analog of vanishing one-loop colored diagram 
of [ 66 ]. In the equations above (ux) = vq f^(nq}o 
is the mean number of magnons per lattice site. Further 
examples, including three-loop corrections to self-energy 
from 7 / 4 “^ and can be found in [ 66 ]. Also, follow¬ 

ing [ 66 ], we omit combinatorial factors in front of dia¬ 
grams and include them directly in corresponding equa¬ 
tions. For example, each of colored contractions from 
equation (15) is accompanied by a combinatorial factor of 


16, which is just the number of permutations of black and 
colored propagators which leave the diagram unchanged. 
We note one, more or less obvious, fact concerning one- 
vertex diagrams: the number of closed loops equals the 
number of ( 71 ^) 0 ’s. This can easily be checked by inspec¬ 
tion of the equations (15)-(17). 

C. Two-loop correction to the magnon self-energy 

To calculate the two-loop correction to the magnon 
self-energy, one must include magnon-magnon interac¬ 
tions described by and Hf\ A direct computa¬ 
tion shows that one-vertex diagrams, as well as two- 
vertex diagrams with both external lines attached to the 
same vertex, simply renormalize the magnon mass [ 86 ] 
(mo —>■ m(T)), while the ’’sunset” diagrams (each vertex 
carrying single external and three internal lines) change 
the geometry of the magnon dispersion. The result is [ 66 ] 




E(fc, icu„) = ^ [A{T) + A{T)B{T)] + 




1 


2 m. 


2mn 25*2 1 2D 


2mn 


'p,<i 




p-q 




Fp.q (i^«) = 


(np)o[l T ( 77^)0 -f {P'k+p—q}o] (^q)o(^fe-t-p—q)o 

ujo{k + p - q) - ujo{p) + uoiq) - iuJn 


(19) 






(^p)o[('^p)o “1“ 1 ] 

f 




where ^{iun) comes from the Matsubara summa¬ 
tion in the sunset diagram, and (7ip)o[(7ip)o + ^]/T 
is induced by two-loop diagrams in which both exter¬ 
nal lines are attached to the same vertex (See [66]). 
Renormalized magnon energies are now easily found 
as W 2 ioop(fc) = wo(fe) - ^W 2 ioop(fe), with (5a;2ioop(fc) = 
lim 5 _j.o ReS(A;, a;o(A;) -|- W). It is seen that magnons, as 
pions in Lorentz-invariant theories [79, 87], remain gap¬ 


less at two loop. In the rest of this section we shall 
exploit perturbation theory to test certain approximate 
treatments of magnon-magnon interactions against two- 
loop result (19). The calculation that led to (19) will 
be referred to as true magnon interaction theory, since 
it preserves spin-rotation symmetry up to and including 
two-loop corrections to the magnon self-energy. 
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D. The Random Phase Approximation 


This well known approximation [12] is usually de¬ 
scribed as the one in which correlations between and 
operators from adjacent sites are neglected. Though 
it is most frequently discussed within TGF formalism, 
it can be easily obtained by replacing S^{t) with site- 
independent average (S^) in the equations of motion 
and commutation relations for The net effect 

of this approximation is the renormalization of magnon 
mass, mo —>■ mo x S/{S^), leading to incorrect descrip¬ 
tion of thermodynamics at low temperatures (the pres¬ 
ence of so-called spurious term in low-temperature 
series for spontaneous magnetization). It should be 
obvious from renormalized magnon spectrum and low- 
temperature series of spontaneous magnetization that 
RPA transforms localized interacting spins into a sys¬ 
tem of lattice bosons, as one can introduce Schrddinger 
field operators {y,t)] = voA{x — y), in terms 

of which 


WRPA(fc) = -z -, 

zmRPA 


S 


mRPA = mo 




( 20 ) 


and the RPA description of 0(3) HFM formally cor¬ 
responds to the free field solution of this bose theory. 
For example, the correlation function {S~ S^) can be 
written as (5'“S'+)rpa = 2{S^)vo{'il}Hx)tlj{x))B.pA = 
2{S^)vo fi^{nk)RPA, where {nk)RPA denotes the Bose dis¬ 
tribution for magnons with energies WRpA(fc). This 
agrees with standard TGF theory [12]. Much less trivial 
is the question of magnon-magnon interactions induced 
by RPA. 

After original papers of Tyablikov and Bogolyubov, a 
number of authors [88-92] succeeded in obtaining the 
same results using different techniques. Where they 
all had failed is to connect RPA results with magnon- 
magnon interactions and internal symmetry of the model, 
since all those works (including the original one of Tyab¬ 
likov and Bogolyubov) are spin-operator oriented. In 
particular, the perturbation theory in [89, 90, 92] is built 
on convenient MFA and fermion representation of spin 
operators, respectively, while [88, 91] develop linked clus¬ 
ter and 1/Z expansion. 

It is only recently [66] suggested that RPA can be ob¬ 
tained in physically more transparent manner, using the 
perturbation theory for lattice magnon fields. The ap¬ 
proximation in question is the replacement of i7eff = 
Hq + -b with Hrpa = Hq — leading to 

the model in which magnon-magnon interactions specific 
to the model with a WZ term are omitted (see [66] for 
details). The one-loop self energy in this approximation 


can be written as 


SRPA(fe) = 



1 ^^0 
S 2mo 


{np)o p-k 


( 21 ) 


where 


RPA 


m 


denotes the corresponding ver¬ 


tex. Of course, this cannot be the final result, since 
(21) implies the gap in magnon spectrum. To get phys¬ 
ically acceptable self-energy, we must impose the con¬ 
strain I := vq fp{np}oP^ = 0. But this is equivalent 
to Vo Jp{np}o'y(p) = vq Jp{np)o = {rix), and the con¬ 
strain actually introduces site-independent mean number 
of magnons, resembling the replacement of the operator 
S^{t) with site-independent average (S^) in TGF theory. 
Renormalized magnon energies become 


WRPA(fc) = Wo(fc) — ERPA(fc) = - - i—(22) 

27710 *8 

which is indeed the low-temperature limit of RPA result. 
If we reverse the steps leading us to the effective Hamil¬ 
tonian of lattice magnon fields, we see that RPA can be 
taught of as arising from the effective Lagrangian 

= I {dtU^u^ - dtU^u^) - 

- ^-n-^ix)ATT'^{x). (23) 

which manifestly violates internal 0(3) symmetry of the 
Heisenberg spin Hamiltonian starting at 0{p^). We 
stress that the reduction of magnon-magnon interac¬ 
tions is is crucial for thermodynamics of 0(3) ferromag- 
nets, for spurious terms in low-temperature series for 
free energy are present whether or not the constrain 
^ofp {np)oP^ = 0 is implemented [66]. 

As noted in the Introduction, RPA is widely used tool 
for studying magnetic systems [12-28, 31-42]. Thus, to 
be able to correctly interpret results obtained from RPA, 
it is of interest to fully grasp the approximation in itself. 
The spin operator-based mathematical machinery of two- 
time temperature Green’s functions makes the physics 
involved less obvious, and this can lead to a wrong con¬ 
clusion. For instance, the authors of a recent paper [93] 
characterize RPA as a theory which ” which does not take 
account of spin correlations or magnon-magnon interac¬ 
tions” . This is at best misleading, since we have shown 
that RPA corresponds to an effective field theory in which 
magnon-magnon interactions arising in the WZ term are 
omitted, but those induced by the unimodular constraint 
are kept. Thus, the difference between LSW and RPA 
comes from magnon-magnon interactions described by 
—H^'^ of (10). The thermodynamic averages in RPA are 

actually calculated using Hrpa = Hq — and not the 
full Heisenberg Hamiltonian, as one may expect from the 
mathematical formalism of TGF theory [12, 94]. 
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FIG. 2. The diagrams constituting magnon self-energy in the Callen approximation. 


E. The Callen Approximation 


1. Magnon Spectrum 


Next we consider the Callen approximation (CA) (See 
[12] and references therein), which may also be viewed 
as a linearization scheme for EOM of spin operators. 
The Callen’s linearization approximately accounts for the 
short-range fluctuations of the order parameter, appar¬ 
ently neglected in RPA. As a result of CA, magnons ac¬ 
quire renormalized mass 


mcA = ™ 


RPA 


1 + 2 


252 


Vo / {nq)cAl{q) , (24) 

Jq 


with unchanged geometry of dispersion relation 

iOCAik) = 


2mcA 


(25) 


{S^)/2S^, usually denoted as a{T), is a phenomenologi¬ 
cal parameter introduced by Callen to get a ’’correction 
to RPA” [95] {a(T) = 0 directly leads to RPA). Since 
the discussion from previous section pointed out on the 
neglection of magnon-magnon interactions generated by 
WZ term as a defining feature of RPA, the question of 
magnon-magnon interactions induced by CA naturally 
arises. Now we examine it in detail. 

First, we observe that the entire effect of Callen’s lin¬ 
earization again reduces to the pure renormalization of 
magnon mass. This means that the magnon-magnon 
interactions are described solely by one-vertex interac¬ 
tions (perturbations), i.e. sunset and other similar dia¬ 
grams are absent as in the case of RPA. Further, using 
(S^) ~ S' — (rix) which is certainly good approximation 
at low temperatures, we see that m-cA contains terms lin¬ 
ear, quadratic and cubic in (riq). According to the gen¬ 
eral structure of one-vertex corrections to the self-energy 
outlined in Sec. Ill, the effective Hamiltonian that ought 
to produce (25) through perturbation theory must con¬ 
tain terms with six magnon operators. However, they 
are absent from (9), since products of six magnon oper¬ 
ators originating from term with spatial derivatives and 
WZ term precisely cancel. This is a special characteristic 
of the effective Lagrangian (8), and terms with six Gold- 
stone fields generally appear in Lorentz-invariant theories 
[51, 79]. The most general effective Hamiltonian with six 
magnon field operator and single lattice derivative is 


He = H^e"^ + (26) 


with 

H'^e^ = AeF'^vo ^ ^ tt{x) ■ 

X 

= BqF'^vo ^ [k'^(x)\^V'^tz'^(x), (27) 

X 

H^e'’ = CqF'^vo ^ 7r^(a;) 7r(a;) • [ 77 ( 0 ;) 7r^(a;)] , 


where Ae,Be and Ce are arbitrary constants. We can 
safely choose Cg = 0, since this kind of interaction is not 
generated neither by daU ■ daU, nor by WZ term of (8). 
Also, we may set Ag = 0 as the the contribution of 
to the self-energy 



(28) 


contains the term oc k'^{nx)'^S which does not appear 
in wcA(fc)- As in earlier sections, denotes the 


corresponding vertex. 

As far as eight-magnon terms are concerned, one eas¬ 
ily verihes that and can not produce desired 
terms in magnon self-energy. The corrections arising 
from are listed in (15)-(17), revealing redundant 
contribution from (17). Also, gives unacceptable 
term oc k^{rix)^S~^. In fact, the only way to reproduce 
vJcA{k) through the perturbation theory is to include the 
new eight-magnon term 


= CeF^vo Y, [^'(*)]' V27r2(a^). (29) 

X 


If the remaining coefficients from (27) and (29) are chosen 
to be Be = 1/32 and Cg = —1/576, the Hamiltonian 

He A = He + + H^'^ + H^^'^ + 77^^=^, (30) 

with suitable constrain on self-energy, yields a;cA(fe) 
through the perturbative corrections. To confirm this, we 
simply evaluate the one-vertex self-energy for diagrams 
depicted at Fig. 2. Since 
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(33) 


(34) 


(35) 


where, as in the case of RPA, I = vq J^(nq)oq'^. The 
absence of external momentum in terms proportional to 
X and makes the self-energy (35) unphysical. These 
arise from two and three loop corrections that don’t in¬ 
clude true interactions generated by WZ term and uni- 
modular constrain thus explicitly breaking 0(3) symme¬ 
try at six and eight magnon field operators (see equations 
(44)-(46) below). Therefore, the contributions to the self 
energy calculated with and must be supple¬ 
mented with additional constrain 1 = 0 for a magnon 
spectrum to display the Goldstone mode. With restric¬ 
tions on X, the self-energy (35) reduces to 


ScA(fe) 


' (Ucc) 

2mo S 


X 

S 


2{nx)X 


I(na,)2' 

53 


,(36) 


where I = vo fq{nq}ojD(q)- The magnon energies in this 


approximation read 


WCA(fe) = -ScA(fc) 

Zttiq 

_ S - (Ux) r S-{nx)q. 
~ 2mo S [ 52 


(37) 


and we recognize (37) as low-energy limit of Callen’s re¬ 
sult (25). From the point of view of an effective theory, 
all terms in magnon self energy arising from symmetry- 
violating parts of the magnon Hamiltonian must be con¬ 
strained with X = 0 to mimic 0(3)-invariant perturba¬ 
tion theory. It is well known that fine tunings of this sort 
are a necessity if the approximation in question is too 
crude to encapsulate fundamental properties of a model 
under study [96]. In the present case it is the replace¬ 
ment of the lattice magnon Hamiltonian (9) with (30) 
that destroys spin-rotational invariance and affects the 






























































10 


magnon self-energy (35). Note, however, that the con¬ 
strain X = 0 should not be enforced on leading term in 
(35), since preserves spin-rotational invari¬ 

ance at four magnon operators. In contrast to RPA, the 
constrain X = 0 merely eliminates the gap and should 
not be interpreted as introducing the site-independent 
number of magnons. 

Thus, Callen’s modification of RPA actually consists 
in adding three more interacting terms to iJupA. The 
one with four magnon fields carries interactions 

generated by the WZ term and preserves spin-rotational 
invariance up to one-loop corrections to the self energy. 
On the other hand, and describe additional 

interactions which are, however, not generated by the 
effective Lagrangian ( 8 ), as seen from (lO)-(ll). The dif¬ 
ferences between ScA(fc), ^RPA{k), and two-loop result 
(19), which are obvious, are thus fully explained in terms 
of magnon-magnon interactions. 

2. Low- Temperature Thermodynamics 


that their contribution to the free energy of 0(3) HFM 
on a simple cubic lattice yields terms whose powers of 
temperature are 5/2, 7/2, 9/2 and so on. Moreover, the 
coefficients of two-magnon terms are the same in CA as in 
true magnon-interaction theory (compare (38) with equa¬ 
tion (71) of [ 66 ]). It is therefore not surprising that CA 
exactly reproduces first three terms in low-temperature 
series for spontaneous magnetization [95], as they are a 
simple consequence of the lattice anisotropies. 

Much more interesting are multiloop diagrams as they 
arise due to magnon-magnon interactions. Out of two- 
loop diagrams, the first one that describes the magnon- 
magnon interaction, and is a potential candidate for 
term, comes from four-magnon parts of Yet, it van¬ 
ishes since contributions from unimodular constraint and 
WZ term precisely cancel as in true magnon interaction 
theory [60, 61, 66 ]. The fact that low-temperature ex¬ 
pansions obtained by Callen [95] do contain term for 
S' = 1/2 should be ascribed to the pure artifact of TGF 
formalism [100]. Thus, the first nonzero contribution to 
the free energy arising solely from magnon-magnon in¬ 
teractions in CA is the Dyson term (proportional to X®), 
the two-loop diagram of [ 66 ]: 


To gain further insight into the CA, we shall con¬ 
sider the low-temperature series for free energy compar¬ 
ing, where necessary, the prediction of CA with those of 
true magnon interaction theory [97]. Starting from the 
magnon Hamiltonian (9) on simple cubic lattice we ob¬ 
tain Hamiltonian density [ 66 , 98] up to p® 

no = ^daTT ■ 5aTT - EfiH (l - Tt'^/2) , 
n^^^ = [tt^tt • Att — 


18 


^ 32 

HW = -hdl-K ■ dl-K 

+ 7 [^l - dlTr^dlir^] 

+ -^ 

16 

= cid^Tv ■ 


(38) 



(39) 


^ „-uHn/T 

T.- - 

n—l 


,5/2 


rj^O 


The next two-loop diagram, including four-magnon ver¬ 
tex of as well as two and four magnon vertices of 
gives X® term (see [ 66 ]). Just as d/ 5 , it occurs both 
in CA and in true magnon interaction theory. Still, this 
is not the leading correction to the Dyson term. 

The next term describing magnon-magnon interactions 
comes from three-loop diagrams. In true magnon interac¬ 
tion theory, this is the three loop diagram of 77 [61, 66 ] 






d^Tv d^TT - — 5a7r®5a7r 

3_ o3_ 

< 5 / 11/2 = -T i' 


(*) 


I{h)T 


11/2 


with external magnetic field along z—axis, formal val¬ 
ues of coupling constants h = X^a^/24, ci = F‘^a^ll2Q 
and A = dada- Given the Hamiltonian density (38) and 
knowing that the loops are suppressed by X) = 3 pow¬ 
ers of momentum (see ( 12 )), we are in the position to 
set up a diagrammatic series for the free energy density 
/ = F/V [ 66 , 99]. The magnon vertices in continuum 
field-theoretic setting will be marked by corresponding 
power of momentum with exception of those generated 
by 77denoted by a dot. 

The diagrams split into two categories [ 66 ]. The first 
one collects all one-loop diagrams with two-magnon ver¬ 
tices which describe lattice anisotropies. It is easy to see 



h = tH/T, (40) 

where, in the absence of external magnetic field 77, 
7(/i = 0) = 5.3367 X 10“® [ 66 ]. As the analysis from 
previous section showed that the CA magnon dispersion 
emerges through the perturbation theory only if interac¬ 
tions described by two-vertex (and higher) diagrams are 
neglected [101], the CA low-temperature series for free 
energy does not contain X^^/^ term of (40). However, 
a term proportional to X^^/^ is present in CA, and it is 
generated by three-loop diagram of 77^^^. Since 
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with first colored contraction appearing four times and 
second one vanishing, we find 


5f 


CA _ 
11/2 ~ 


2 

/S' 

.9/2 ^9/2 


£2 


1 {27rf 

2^ 7 ^ 3/2 

.n—1 


X 


n—1 


-fiHnIT 


,5/2 


2111/2 


(42) 


By comparing (40) and (42) for H = 0, we that CA over¬ 
estimates the leading correction to Dyson’s term nearly 
by a factor of 39. Obtaining expression (42) within TGF 
formalism is somewhat tedious, but rather straightfor¬ 
ward within magnon perturbation theory, i.e. when ap¬ 
propriate set of magnon-magnon interactions have been 
identified. 

Finally, we note that Hamiltonian (38) can be obtained 
from the Lagrangian 

C = £[2] + £[4] + £[6] (43) 


where 

£'^1 = ■ d^U + 


+ - ^daTV^da'n'^ 

32 576 

- (^^Tt) • 5„7r, 


(44) 


^1 




- ^ [dl • dl-K - dlir^dlTZ^] , 


(45) 


£[®1 = c^U • dldlU - ^ 

Cl 


1 

18' 






(46) 


It is seen that each £[^], £W and £[®1 contain terms which 
manifestly violate internal symmetry of the Heisenberg 
model, just as in the case of RPA. 

The explicit expressions from Sec. HI D and Sec. HI E 
should supplement imprecise statements relating RPA 
to CA found in the literature and dating as back as 
Callen’s original paper [95]. For example, Tahir-Kheli 
introduces CA as an approximation which takes ’’into 
account the fluctuations of around its average (S'^)” 
[102]. Further, it is claimed in a recent paper [93] 
that ”CA takes, to some extent, account of magnon- 
magnon interactions”, without actually specifying their 
type. Also, the point of view in standard and up-to-date 
reference on TGF method [12] is that CA ’’takes some 
higher-order correlations into account”. By striping the 
issue of RPA/CA relationship down to the problem of 
magnon-magnon interactions, we have shown that effec¬ 
tive Hamiltonians (i.e. Lagrangians) that generate RPA 
and CA through perturbation theory violate 0(3) sym¬ 
metry of the Heisenberg model at 0{p^). A difference is, 
of course, that manifest symmetry breaking in CA takes 
place at six and eight magnon terms developing devia¬ 
tions from true magnon interaction theory at three-loop 
corrections to the free energy, while RPA breaks down 
already at two loop. 


F. Comparison with quantum Monte Carlo 
simulation 


We conclude the analysis in this section by compar¬ 
ing predictions of effective field theories for RPA and CA 
with quantum Monte Carlo (QMC) simulations and true 
magnon interaction theory for the 0(3) ferromagnet free 
energy. This will allow us to infer the influence of RPA 
and CA-type magnon-magnon interactions on thermody¬ 
namics beyond leading order expansions in T discussed 
so far. The usefulness of the lattice EFT pursued here 
will become clear when we look at results away from 
low-temperature sector, i.e. when the lattice structure 
becomes resolved by magnons. While equations given 
bellow apply equally for arbitrary lattice, exchange in¬ 
tegral J and localized spin S, all numerical calculations 
are conducted for the S' = 1/2 Heisenberg ferromagnet 
on a simple cubic lattice with J = 1. The parameters 
S = 1/2 and J = 1, together with the lattice type, define 
a concrete realization of the Heisenberg ferromagnet. On 
the other hand, the parameters of EFT are F and E and 
are related to those of original Heisenberg Hamiltonian 
by matching [see the discussion bellow equation ( 11 )]. 

The free energy (per lattice site) corresponding to the 
non-interacting theory is calculated with FIq of (9). This 
is basically a LSW result and is given by 


Flsw 

N 





4 _ Q-‘-^o{k)/T 


(47) 


To make comparison with QMC easier, we have normal¬ 
ized the non-interacting free energy by adding Fq/N = 
—3/4, the ground state energy of S = 1/2 and J = 1 
Heisenberg ferromagnet on the simple cubic lattice. We 
note that Flsw includes all one-loop diagrams (in the 
sense of derivative expansion) permitted by the lattice 
symmetry. 

The two-loop correction to Flsw in true magnon in¬ 
teraction theory (TMIT) is found in [ 66 ] to be 


N 


1 ^^ 

S 2mo 4D 


Vo 


/ {np)o P ‘ 
Jv 


1 2 


(48) 


and the two-loop correction corresponding to RPA is 



with RPA vertex defined bellow equation (21). According 
to the discussion from Sec. HIE, CA type of magnon- 
magnon interactions induce two, three and four-loop sin¬ 
gle vertex contributions to the ferromagnet free energy. 
While the two-loop term is the same as in TMIT [equa¬ 
tion (48)], the three and four loop corrections are given 
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r[K] 


FIG. 3. The influence of magnon-magnon interactions on the 
free energy of S' = 1/2 and J = 1 Heisenberg ferromagnet on 
the simple cubic lattice. The purple (top) curve represents 
the non-interacting model (LSW), while blue and black one 
represent TMIT and QMC results. The orange and green 
lines correspond to the models with interactions of CA and 
RPA type. The value of critical temperature Tc = 0.839(1) 
K is taken from [105]. 


by 



_ 1 2 

2S^ 2mo ° 



{nq)oP-q ■ (51) 


Finally, the simulation was performed using stochas¬ 
tic series expansion within quantum Wang-Landau algo¬ 
rithm [103] based on ALPS libraries [104]. The cutoff is 
set to 5 X 10^, which, for the SC lattice with 10^ sites, 
yields reliable results if T > 0.05J (with J = 1). By let¬ 
ting the number of Wang-Landau refinement steps to be 
18, the total number of sweeps becomes ~ 10^° and the 
error bars are much smaller than the line used to display 
the simulation results. 

Results of calculations are collected on Fig. 3. First, 
it is obvious that all lattice magnon models (including 
LSW) develop only small deviations from QMC almost 
up to T = 0.5 K. This clearly confirms our remarks 
about weakness of magnon-magnon interactions even in 
the case of S' = 1/2 since the critical temperature is, 
from QMC simulations, known to be Tq ~ 0.84 K [105]. 


It is also evident that TMIT produces results closest to 
QMC in wide temperature range, as expected, demon¬ 
strating efficiency of the lattice magnon theory. Proba¬ 
bly the most important observation concerning practical 
applications and interpretations of RPA and CA is that 
both approximations are based on interacting theories 
which, thanking the loop structure of corrections and 
corresponding Hamiltonians, overestimate the influence 
of magnon-magnon interactions [compare with equation 
(42)]. This is the main reason why self-consistent de¬ 
scriptions of real compounds based on RPA/CA, as a 
rule, predict higher values of the exchange integral J 
than LSW or non-linear spin-wave theories based on bo¬ 
son representations [106]. For example, the non-linear 
spin wave theory in combination with experimental data 
for 3D ferromagnet CrBra predicts the intra- and inter¬ 
layer nearest neighbor exchange integrals J = 8.25 K and 
J' = 0.497 K [107]. In contrast, self-consistently deter¬ 
mined exchange integrals within RPA are J = 12.38 K 
and J' = 1.0 K [8]. Of course, the strength of magnon- 
magnon interactions in RPA and CA at higher temper¬ 
atures is masked by the special form of self-consistent 
equations for correlation functions. Namely, RPA and 
CA solutions in TGF formalism share many features with 
simple mean-field description at higher temperatures, ex¬ 
emplified by the expressions for the Curie temperature 
and the value of critical index /3. As a result, both RPA 
and CA predict somewhat higher value for the critical 
temperature ps 0.989 K and m 1.326 K for 

J = 1 and 5 = 1/2 ferromagnet on a simple cubic lat¬ 
tice). Finally, we note that the similar conclusions con¬ 
cerning low-temperature behavior of 0(3) ferromagnet at 
low temperatures were reached in [108]. 


IV. PERTURBATION THEORY FOR TYPE A 
MAGNONS 


Now we turn to the perturbation theory for type A 
Goldstone bosons. Basic diagrammatic rules remain the 
same, but the definition of thermal propagator changes. 
Also, we shall need to include more than one lattice 
Laplacian per vertex as well as couplings beyond near¬ 
est neighbors. All this can be achieved within formalism 
of colored propagators described in Sec. IIIB. 


A. Kondo-Yamaji Equations 

The approximation of Kondo and Yamaji (KYA) rep¬ 
resents an extension of the EOM method in the sense the 
linearizations are conducted in equations containing S^. 
The long-range order parameter {S^) does not explic¬ 
itly enter the magnon dispersion, making this approach 
suitable for low-dimensional systems [43, 44]. It is also a 
starting point in further developments of TGF formalism 
[19-21, 44]. Here, the focus will be on one-dimensional 
ferromagnet. Strictly speaking, effective field theory is 
not directly applicable to one-dimensional 0(3) ferro¬ 
magnet in the absence of external magnetic field, since 
magnons acquire nonperturbatively generated gap. How- 
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ever, a careful analysis from [64, 65, 109, 110] reveals 
that first few terms in low-temperature series for free en¬ 
ergy are well defined even in case of zero external field, 
and KYA yields results in agreement with thermal Bethe- 
ansatz in this extreme low-T sector [ 111 ]. 

The Kondo-Yamaji equations (KYE) for magnon dis¬ 
persion 

^KY{k) = J {[1 — cosfca][l — Cl -|-C 2 — 2 ci cosfca]}^/^, 
_ r noth 

Cn = Wvo J 2a;KY(p) ~ 


with IT ~ 2 J and cq = 3/2 at low temperatures, 'y{p) and 
given in (2) and (13) for D = 1, are derived from the 
assumption of vanishing long range order (LRO) [43]. It 
turns out that is extremely useful to exploit the constrain 
Co = 3/2 directly and to write KYA dispersion as 


ujiy{k) = [3ho-ci] 


; 2 2 
k C 


a^Wvo 


coth 


t^KY(p) 

2T 


2 wky(p) 


(53) 

[p^]^ + fc^fc^a^c^ci. 


c 2 = 




We shall demonstrate bellow that KYE (52) i.e. (53), can 
be obtained from the effective Lagrangian description of 
a ferromagnet chain if E = 0 is set in ( 8 ). When E = 0 
the effective Lagrangian must include a term with two 
temporal derivatives and it takes pseudo-Lorentzian form 


CeS = ^dtU^dtW - (54) 


where, as before, U denotes the unit vector field and 
and 7 ^ are constants satisfying c^ = jF"^. In con¬ 

temporary nomenclature, this Lagrangian describes the 
dynamics of type A Goldstone bosons [57, 69]. 


B. Hamiltonian and perturbation theory for NN 
lattice model 


To develop the perturbation theory that yields KYE, 
we construct the lattice magnon Hamiltonian. Since the 
Lagrangian (54) contains [9tC/]^, the Hamiltonian can¬ 
not be obtained as in preceding sections. Instead, we 
follow [112] to get the interaction picture Hamiltonian. 
By putting it on a spatial lattice (chain), we find the free 
part 



with (thermal) propagator 

-y^Tx- Ty) = (T {'K°‘{x,Tx)TT^{y,Ty)])o (56) 


^ab 

T 


E 


n— — oo 


Wo(9)+w2 


free-magnon dispersion 


(57) 


and nearest-neighbour (NN) interaction piece 


Hint 



Et^^V^tt^ - ^vo 


^[7r-5^7r]^ (58) 


where 7r“ denotes the interaction picture magnon field 
and n“ is the corresponding canonical momentum. 

To calculate one-loop self energy for interaction (58) 
we use diagrammatic rules outlined in Sec. HIB. While 
the action of on propagator lines merely introduces 
a factor of g ^ for external or loop momenta, propaga¬ 
tors with temporal derivatives must be handled with care 
[112]. We find 

Dlb = (T {dr^TT‘^{x,Tx)TT\y,Ty)})o 



D-y ijli ^n) 


^\q{x-y)-\uJr^{Ta^-Ty) 


D 


— T 

ah 


{F {-K'^{x,Tx)dryTr^{y,Ty)])o (59) 

^ n=-oo-^q 


D 


TT 

ah 


(T {dr^'^‘'[x,Tx)dryTT^{y,Ty)})o 
P n=-ooJq 


where 


Dr{q,UJn) = -D_r[q,^n) 


Hrrhli a^n) 


-^l{q) 

Wo(g) +w 2 


-i^n 

Wo(g) +w 2 


(60) 


Using propagators listed in (59) together with (56) and 
action of lattice Laplacians, we obtain one-loop self en¬ 
ergy 


E(fc,a;) = 


coth 


2ujo{p) 


= 7^vo 


coth 


t^aip) r„2 


2T 


Jp ‘2uJoip) 

and thus 

u^{k) = uj^{k)-i:{k) 


—k^p^ + ^ 


= k^c^ 


1 - F^vo 


coth 


■^oIp) „2 

2T “ ; 


2 wo(p) 2 


(61) 


(62) 


where E(A:) = E(fc, a;o(fc)) is the on-shell self energy. If we 
set = 3IU ~ 6 J and interpret (62) as a self-consistent 
equation for w(fc). 


w^(fc) = (?k^ 


1 - 


= [I - 



3co], 


coth ^2 

2uj{p) 2 ^ 


(63) 


we see that the lattice model (58) with NN coupling ac¬ 
counts only for the leading term in EKY(fc) listed in (53). 
Obviously, the full KYA self energy can be obtained only 
with additional interactions in the one-loop approxima¬ 
tion. 


wo(g) = c^/^ = cq, 
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C. Perturbation theory for KYA model 


1. Colored diagrams for NNN interaction 

Two additional conventions will enable us to imbed 
next-nearest neighbor (NNN) interactions, as well as in¬ 
teractions containing more than one lattice Laplacian 
into the framework of colored diagrams. They will also 
make diagrammatic rules suitable for potential higher- 
loop calculations. First, propagators affected by V^ 2 ) 
be denoted by a double colored line. Second, the presence 
of V^ 2 ) 01 ' be represented by colored edges of the 

pentagon, which denotes a vertex. For example, typical 
interaction containing both V^ 2 ) NN Laplacian is 
V( 2 ) 7 r)(V^ 7 r • V^tt). Four distinct colored 
contractions contribute to one-loop self energy generated 
by V. These are 



since each of diagrams listed above appears four times. 


2. KYA self-energy 


As a first step, we introduce the second-neighbor cou¬ 
plings in (55) and (58) by the replacement —>■ 

, where 


:= = 



(65) 


V^ 2 ) being defined in (5) and the strength of the next- 
nearest neighbor (NNN) coupling G will be determined in 
perturbative expansion. Much of the results obtained for 
the NN model still hold. In particular, the propagator 
and the one-loop self energy retain their general form 
given in (56) and the first line of (61), respectively, with 
free dispersion being 


nl{k) = 


k'^ -f 


c^=j^jF^. ( 66 ) 


Further, let us now introduce additional interaction 
terms 


= A(2)c2uo ^(tt • VV)(V27r • V^tt) 

X 

^(vV • ■ V^tt) 

X 

H^^I=Ac^voJ27t-V^7t (67) 

X 


where and serve as counter-terms, i.e. they 
absorb renormalizations of and couplings. 

Individual contributions of (67) to the self-energy are 
evaluated using colored diagrams, as outlined above. To 
remain within KYA, we discard a term proportional to 
k^k^k^ and arrange the remaining ones in two groups: 
the ones proportional to k"^ and those proportional to 
The corresponding self-consistent equation for 
magnon energy reduces to KYA form (53) if A = 251/48, 
B + 0"^ = 7^(2) = a®lT/16, = -la^WjM, 

= c?W and 7 ^ = c^F^. From the compact form 
of magnon energies [equation (52)], one readily shows 
that Cl = C 2 = 1/2 at T = 0 K so that the magnon 
dispersion becomes identical with standard LS’W result, 
ui{k) = J{a^/2)k^. This is to be compared with fig oc k. 
Thus, in contrast to the perturbation theories for type 
B magnons from preceding section, interactions substan¬ 
tially change the magnon dispersion law in the present 
case. The possibility that interactions between Goldstone 
bosons may have such a strong impact on the geometry 
of their dispersion was anticipated only recently [113]. 
To the best of author’s knowledge, calculations from this 
section offer the first explicit realization of such a mech¬ 
anism. At finite temperatures, we find 


Cl 


I3/2 — 


- - - 2T^, 

2 2 

C(3/2) 


C2 


973/2 


leading to 


( 68 ) 


\k) ^ U'^a^T^k^ + 


2„4 r 


J^CL 


1 - 373/2r = 


k^ (69) 


in agreement with [43]. 
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T[K] 

FIG. 4. The influence of interactions on the free energy 
of S' = 1/2 and J = 1 Heisenberg ferromagnet on the 
chain lattice. The purple curve represents the standard non¬ 
interacting model (LSW), while the black one denotes QMC 
results. The orange line corresponds to the model with inter¬ 
actions of KYA type. 


The results shown at Fig. 4 compare KYA [equation 
(70)] to LSW and QMC. We see that, in contrast to 
naive LSW, KYA result agrees with QMC rather well in 
wide temperature range. Further, the present exposition 
reveals the mechanism behind KYA; similarly to RPA, 
magnon-magnon interactions generated by WZ term do 
not contribute to KYA. Instead, the dynamics of a fer¬ 
romagnetic chain is governed by NN and NNN couplings 
for type A magnons, combined in a such manner that 
one-loop corrections yield spectrum usually encountered 
at type B systems at T = 0. Even though the KYA devi¬ 
ates from QMC only slightly, which makes it a relatively 
good method for qualitative description of a S' = 1/2 fer¬ 
romagnet chain, it introduces spurious terms arising from 
magnon-magnon interactions in the low-temperature se¬ 
ries for free energy already at leading order, i.e. at 
[65, 115]. 


V. SUMARRY 


3. KYA Thermodynamics and Quantum Monte Carlo 
Simulation 


Now we discuss the low-temperature series for free en¬ 
ergy in self-consistent KYA. According to (69) and (54), 
magnon fields admit standard type A Fourier decompo¬ 
sition including both creation and annihilation operators 
which produces so-called zero-point term in expression 
for free energy [99]. The zero-point term is temperature 
dependent in the present case [see (52)] and, since there 
are two type A magnons, we find 


N 


El 

N 


Tvq 



+ Vo 


wky(p), 


X _ q-‘^o{p)/T 


(70) 


where Fq/N = — 5/4 provides a proper normalization at 
T = OK. To see that (53) indeed reproduces all main fea¬ 
tures of low-temperature series for free energy in KYA, 
we calculate (70) with the full Kondo-Yamaji dispersion 
[43] and the approximate one [equation (53) and solu¬ 
tions (68)]. The first two terms of these series match 
completely. While the first one, proportional to 
is rather expected [114], the second one is proportional 
to and describes the effects of interactions. The 
term is twice the one obtained by thermodynamic Bethe 
ansatz equations [115] since KYA describes 0(3) ferro¬ 
magnet in terms of type-A magnons (a similar problem 
with term arises in SBMFT [7]). The direct applica¬ 
tion of TGF method yields (7/3)T^ [111] instead of 2T^ 
and this discrepancy may be classified as another artifact 
of TGF formalism. 

To check the validity of KYA result for free energy of a 
S' = 1/2 ferromagnetic chain at arbitrary temperatures, 
we test it against QMG simulation. The simulation was 
conducted using Wang-Landau algorithm and ALPS li¬ 
braries [104]. Results reliable up to T « O.OIK are ob¬ 
tained be setting the cutoff at 3 x 10"^. 


Linearization of equations of motion (EOM) for spin 
operators is a popular and important tool for studying 
magnetic systems. These include not only insulating ma¬ 
terials with localized spins but also itinerant electron 
systems where Heisenberg Hamiltonian arises through 
a mapping on an effective model. If the parameters of 
linearization are chosen carefully, i.e. if an appropriate 
’’decoupling scheme” is used, this method will provide 
reliable results in many cases. Eor instance, solutions ob¬ 
tained with standard linearizations obey Mermin-Wagner 
theorem and agree with Monte Garlo simulations and ex¬ 
perimental results. However, in choosing parameters of 
lineaarization one must rely on physical intuition since 
approximations are not controlled perturbativley (equa¬ 
tions of motion contain the Heisenberg-picture spin oper¬ 
ators) nor do magnon operators appear at any stage, so 
all linearizations are generally considered to be ad hoc. 
This raises several questions concerning proper interpre¬ 
tation of results obtained by EOM. For example, it is 
hard to identify types of interactions generated by vari¬ 
ous linearizations and their influence on thermodynamic 
properties, or to separate interaction-induced effects from 
those related to the choice of dynamical degrees of free¬ 
dom. Because complete understanding of EOM method 
is a prerequisite for its successful application, these issues 
need to be resolved. 

On the other hand, effective field theory (EET) is a 
powerful tool for handling dynamics of Goldstone parti¬ 
cles as it provides systematic (and model independent) 
method for describing interaction effects. From this point 
of view, 0(3) ferromagnet is a rather interesting sys¬ 
tem since magnon-magnon interactions come in two cat¬ 
egories: the ones induced by the unimodular constrain 
and those generated by the Wess-Zumino (WZ) term. 
The main advantage of EFT over EOM is that magnon- 
magnon interaction terms appear explicitly in the Hamil¬ 
tonian (Lagrangian) making the dynamics of magnons 
completely unrelated to the local su(2) algebra of spin 
operators that define the Heisenberg Hamiltonian. 

The present paper deals with interaction effects in- 
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duced by the three most commonly used linearizations: 
random phase approximation (RPA), Callen approxima¬ 
tion (CA) and Kondo-Yamaji approximation (KYA). It 
is shown that all unique properties of linearizations men¬ 
tioned above originate in different treatments of WZ 
term, where type B (A) magnons are described by the 
model with (without) WZ term. Since all linearizations 
discussed in the present paper respect discrete structure 
of the Heisenberg Hamiltonian, we employ lattice regu¬ 
larization that preserves these symmetries. Perturbative 
calculations are conducted with the aid of colored di¬ 
agrams which are particularly suited for theories with 
derivative couplings. 

As the simplest of type B models, RPA describes the 
0(3) ferromagnet with one-loop interaction theory in 
such way that interactions specihc to the WZ term are 
neglected [see Sec. HID]. This simplification induces vi¬ 
olation of 0(3) symmetry at leading interaction terms of 
order and eventually produces spurious terms in low- 
temperature series for free energy. Even though the the 
0(3) symmetry is explicitly broken by interaction terms 
in effective Hamiltonian, the constrain vq JpP^{np)o = 0 
in the one-loop self energy eliminates the gap at the 
same time introducing site independent average number 
of magnons. Thus, opposed to some claims found in the 
literature [see the discussion in the last paragraph of Sec. 
HI D] , magnon-magnon interactions are accounted for in 
RPA. These are interactions induced by the unimodu- 
lar constrain, i.e. they are generated by the geometry of 
coset space 0(3)/0(2) = 5^. 

The next approximation considered in the text is CA, 
viewed by many as an improvement over RPA. Unfortu¬ 
nately, the claims of superiority of CA over RPA have 
been supported only by some vague statements concern¬ 
ing additional correlations or interactions induced by this 
linearization scheme [see the discussion at the ending 
of Sec. HIE2], or by a subsequent analysis of low- 
temperature series for spontaneous magnetization. How¬ 
ever, a posteriori analysis of low-temperature series can 
not distinguish between interaction effects and pure ar¬ 
tifacts of the formalism and therefore it can not be con¬ 
sidered as a complete one. In contrast, analysis based 
on the effective field theory presented in Sec. HIE re¬ 
veals that, regarding magnon-magnon interactions, CA 
does contain certain improvements over RPA. Concretely, 
four magnon interactions arising from WZ term appear 
in CA through one-loop corrections in the magnon self¬ 
energy [see Sec. IIIElj. However, additional spurious 
six and eight magnon interaction terms included in CA 
violate 0(3) symmetry at through two and three loop 
corrections to the self energy. Therefore, the constrain 
Vo jpP^{np)o = 0 must be imposed on those contribu¬ 
tions to the self energy that arise from symmetry vio¬ 
lating terms, just as it is in RPA. Also, one should note 
that all corrections to the self energy in CA come from 
one-vertex diagrams. Because of that, magnon energies 


in CA, as well as in RPA, posses the same geometry as 
LSW solution (that is, uj{k) oc fc^). Alternatively, as de¬ 
scribed in Sec. HIE2, magnon-magnon interactions in¬ 
duced by WZ term and unimodular constrain constitute 
two-loop corrections to the free energy in CA which gives 
correct results at order T® (this is the leading term in 
low-temperature series generated by interactions). How¬ 
ever, symmetry violating terms produce deviation from 
true magnon interaction theory at order (that is, 

at three-loop corrections to the free energy). 

In addition, the applicability of lattice field theo¬ 
ries for type B magnons in wide temperature range is 
demonstrated by comparison with quantum Monte Carlo 
(QMC) simulation [see Sec. HIE]. Excellent agreement 
between QMC and LSW theory combined with one-loop 
correction to the free energy [equations (47) and (48)] 
additionally supports lattice EFT as a right method for 
interpretation of EOM results. 

Further, we discuss KYA and show that its physics 
near T = OK can be reduced to the self-consistent per¬ 
turbation theory for type A magnons. Therefore, WZ 
term and magnon-magnon interactions that it generates 
do not enter KYA, making this approximation substainly 
different from RPA and CA. Still, at T = OK, KYA re¬ 
produces standard dispersion for ferromagnetic magnons, 
uj{k) oc k^. This property of KYA is usually explained 
with phenomenological vertex parameters which account 
spin-spin correlations and maintain rotational invariance. 
In Sec. IV we show that these characterizations of KYA 
may be replaced with self-consistent one-loop perturba¬ 
tion theory for a model of type A magnons with explicit 
identification of interaction terms. In particular, KYA 
Hamiltonian contains lattice Laplacians that connect first 
as well as second neighbors on a chain lattice and these 
additional interactions play a key role in producing char¬ 
acteristic dispersion at T = OK. Analysis of the free en¬ 
ergy reveals that KYA low-temperature series for free 
energy deviates from rigorous results already at leading 
term which is a consequence of interactions. However, 
at the same time, KYA result differs from QMC only 
slightly in wide temperature range [section IV C 3] and 
provides qualitative description of the low-temperature 
thermodynamics. 

A detailed analysis based on lattice magnon fields pre¬ 
sented in the paper revealed physics behind RPA, CA 
and KYA in case of 0(3) ferromagnets. All three approx¬ 
imations are reduced to equivalent systems of interacting 
magnons thus providing better understanding of widely 
used theoretical tools. 
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